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ABSTRACT 



Adaptive control system techniques have been utilized to inves 
tigate identification of the dynamic equations of operating control 
systemso The identification has been achieved by using model modu 
lation techniques to drive dynamic models into correspondence with 
operating control systems » The system identification then proceeded 
from examination of the model and the adaptive loop. 

The model modulation techniques applied to adaptive control 
systems are breifly discussed These techniques applied to simulation 
studies on an analog computer and a digital computer, and applied 
directly to a position servo-^mechanism , are then discussed The 
discussion presents data showing the accuracy of identification of 
the dominant characteristics of the dynamic equations when the order 
of the equation is either known or approximated It was found that 
these dominant roots could be found within an average accuracy of 57o 
for both complex phase angle and undamped natural frequency 

The authors wish to express their appreciation to Professor 
Richard Co Dorf of the U- S. Naval Postgraduate School, for the 
suggestion of the model modulation method and for his assistance 
during this investigation. 
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CONTROL SYSTEM IDENTIFICATION THROUGH MODEL MODULATION METHODS 



Chapter I 
INTRODUCTION 



Genera I 

The last three decades have seen control system technology grow 
from a series of rules of application to a full regime encompassing 
most of the disciplines of engineering. The study of control systems 
has progressed deeply into analysis and synthesis of linear and non* 
linear control systems. There are, additionally, offshoots of basic 
systems with most exhaustive methodology: examples are the sampled 
data, hybrid analog-digital and multi-loop self adaptive systems. 

The problem of control system identification led to some of 
the original investigation in the control field. Despite the increase 
of knowledge in the field of control engineering, system identification 
has remained a major problem. 

The intent of this project was to investigate a method of control 
system identification utilizing techniques from the growing technology 
of the self-adaptive control system. One restriction was laid upon 
the project at its inception. This was that the method of system 
identification remain independent of the system to be identified 
Ideally then, an operating system should be identified solely through 
examination of the input and the response. To this end a control 
system identifier has been modeled and tested on an analog computer, 
through digital computer analysis, and upon a position control system. 



Control System Identification 

The characteristics and performances of the fundamental components 
of feedback control systems are generally amenable to systematic analysis 
The control system, when composed of these fundamental components, may 
be easily synthesized and analyzed. Components, however, when removed 
from the theoretical state do not exhibit the invariant characteristics 
attributed to them. Furthermore, control systems are often synthesized 
on the basis of transfer function manipulation. This synthesis occurs 
under the assumption of negligible loading of one component system by 
another. This assumption is too often fallacious. In either or both 
of these situations the control system engineer is confronted with a 
control system which does not conform to previous analysis. And if 
this analyzed control system is a portion of some larger process, it 
may be necessary to reevaluate the control characteristics and perform- 
ance, in other words to identify the system. 



In some extremely variant systems the processes of system identi- 
fication are utilized to serve as controllers in a self-adaptive control 
system. The most notable of these cases are the self-adaptive control 
systems currently in use and being further developed for aircraft 
installation. The NASA X-15 is an example. In this installation the 
response characteristics are examined and the results are used to vary 
flight control loop gain to maintain near invariant control response. 

The equipment used for this process is quite complex, but is justified 
by the operating environment of the vehicle. 
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For less complex control systems or in less complex environments, 
the problem of control can be solved by more conventional means, without 
recourse to dynamic self-adaptive system technology. These less complex 
systems require accurate identification of the process dynamics. Once 
the transfer function of a process has been evaluated this evaluation 
can serve as a basis for design of the remainder of a system or for 
readjustment of the process through compensation techniques. 

Identification Methods 

Truxal writing in reference (1) has placed the methods of control 
system identification into four broad categories: 

a. Parameter identification. 

b. Frequency response methods. 

c. Impulse response methods. 

d. Direct evaluation of system differential equations 
of motion. 

Some of the difficulties in one-by-one determination of model 
parameters have been mentioned. Additionally, this method is possible 
only when the detailed operation of a process is reducible to funda- 
mental differential equations in which the primary parameters may be 
included and the negligible ones discarded. This ideal situation 
rarely exists for complex control of complex processes. 

As frequency response techniques are well-known, a description 
of them has been omitted. These procedures are extremely important 
in general control design. This is the most common method of iden- 
tifying the transmittive dynamics of a control process. Many working 
control systems however, cannot be subjected to variable frequency 
oscillations for reasons of system complexity, operating economy, or 
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because of the extensive time and effort required for determination 
of the gain and phase characteristics. Two additional problems occur. 
First, in some circumstances, frequency response testing should not 
be superimposed on regular operation, because testing would interfere 
with designed operation. Secondly, a working system is subjected not 
only to its driving input, but also to noise occuring during the op* 
erating process. The extraneous disturbances of normal operation 
would make the interpretation of frequency test results difficult. 



The impulse response can serve to completely identify any linear 
time-invariant two port system. The various impulse response techniques 
are utilized in many of the more sophisticated adaptive control systems. 
Truxal has cited three primary techniques for measuring impulse response. 
The first and most direct method is to excite a system with an impulse 
and measure the response. Second, the system can be excited by a 
known signal. For an initially inert system, the impulse response 
g^^^ can be determined from the input, output, through 

the convolution integral 

11 

Finally the impulse response can be measured by application of a small 
random signal to the process. If white noise is used as the input, 
the cross correlation function between the output and the noise is the 
impulse response. 



‘^(i) / J(r) -Tl 



The attempt to utilize the impulse response in identification is 
complicated by three problems. Noise has a much more delerterious 
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effect on the output than in the case of frequency response tests. Non- 
linearities are much more difficult to detect than in frequency response 
tests. Finally the mechanization necessary to determine the impulse 
response is almost prohibitive for use in system identification. When 
the impulse response method is used in adaptive system control the 
engineer is generally interested not in the full identification, but 
in identifying only some characteristic of the response. Quantitative 
measurements are made to serve as sufficient identification to allow 
adaptive correction. As an example, the number of oscillations during 
a specific period of a system response would identify the impulse 
response. These approximations are not sufficient for full identifi- 
cation, however. 



The evaluation of the differential equation for a linear process 
offers an additional method of identification. The differential equation 
of a process might be expressed as an equation of input output 

c^^^ in the following form: 



oco + (i.ic + 0, = + 00 0 +. b.r, 

at at ot .c/1 



The a^ term may be evaluated from the nature of the zero frequency or 
steady state gain. If the order of the differential equation is known, 
and the nature of the right hand side is known, the other ^ co- 
efficients may be evaluated by generalizations such as minimizations 
of square of the error e, where e is defined: 



e 




Tf'* 



4 



0 o 0 4 






a c 




1.3 
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The a,, coefficients are evaluated through use of this error criterion 

N“ 1 

by generation of successive differentiations of the response. Unfor^ 
tunately, searching over an (N-1) space for a minimum is a problem of 
potential difficulty. General problems also arise in settling at 
local minimums and in error system stability. 

Many of the problems of systematic determination of the differ- 
ential equation coefficients are alleviated by the use of system-model 
reference. This is a natural variant to the adaptive control system 
technology. The exploitation of this method involves utilization of 
the difference between system and model response to the forcing input. 

This difference has been handled by several methods in present technology. 
In general self adaptive processes, the difference as operated upon 
by an adaptive technique is used to drive the system into operating 
correspondence with the invariant response of the model. 

This investigation was based on a variation of this method. During 
this investigation adaptive techniques have been used to drive the 
model into operating correspondence with the operating system. The 
characteristics of the model differential equation have been modified 
by the adaptive loop. These characteristics then have been used to 
accurately define the coefficients of the system differential equation. 
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Chapter II 



MODEL REFERENCE ADAPTIVE SYSTEMS 



Introduction 

A major advantage of the model reference adaptive system is that 
direct measurement of the operating system differential equation 
parameters is not required. This factor results in considerable 
economy in mechanization of the adaptive components. As previously 
mentioned there are in existence various methods of self adaptive 
control. However, in many of these systems the adaptive loop mech- 
anization requires the assumption of the second order nature of the 
system. An example of this is Osder's system detailed in reference (1). 
For this system, the sign reversals in the output oscillations are 
counted and the count used to correct the second order zeta. Most of 
the impulse response self adaptive systems also require assumption of 
a basically second order system to keep the equipment within allowable 
limits. The model reference system relieves this restriction. The 
model reference adaptive system can be utilized for higher order systems 
about which the second order assumption is invalid. Additionally, 
current literature presents several examples in which model reference 
allows the control of several parameters within the operating system. 

References (1) and (2) present several methods of model reference 
adaptive control systems which have been proposed or developed. The 
variants of these systems include the Ham-Lang conditional control 
system, which in analysis is passive adaptation through model reference. 
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The other extreme is perhaps, Braun's method, which utilizes an ortho- 
normal spectrum analyzer to solve the identification problem. 

Parameter Perturbation 

Perhaps the most straightforward adaptive mechanization results 
from the use of model reference with perturbation techniques. McGrath, 
Rajaraman and Rideout in reference (3) have proposed and examined 
various parameter perturbation model reference schemes,. Figure 1, 
taken from their discussion presents two variants of the parameter 
perturbation scheme . 



Fig. 1 

PARAMETER PERTURBATION METHODS 








The adaptation processes in these systems are driven by the error, e, 
between the model and system. This error is used to form a performance 
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criterion, which for this example was integral of error squared, J e d t. 
This criteria contains a component of the frequency at which a parameter 
of the System In Fig. 1(a) is being perturbed. The amplitude and 
phase of this component give the magnitude and sign of the signal 
obtained through perturbation and Integration. This signal is then 
fed back negatively to reduce the error. The adaptive loop mechan* 
ization is illustrated in Fig. 2. 



Fig. 2 

PARAMETER PERTURBATION ADAPTIVE LOOP 




For analog simulation electronic multipliers have been utilized for 
parameter correction insertion. In operating systems servo'-mul ti- 
pliers could be utilized. 

The system of Fig. 1(a) suffers from the disadvantage that the 
perturbation appears in the output. The system of Fig. 1(b) avoids 
this problem through model parameter perturbation. This second method, 
however, requires that the model be nearly the same form as the system. 
Both methods, moreover, require phase correcting networks to correct 
the phase shift between parameter perturbation signal in the model and 
the oscillator signal. Through utilization of further model forms 
with parameter perturbation rhe authors demonstrate that it is possible 
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to prepare a configuration which is both system and signal adaptive, 
using multiple models with perturbation of either model or system 
parameters. Rajaraman in a later article, Ref. 4, has shown the 
detailed operation of the signal and system adaptive process. In 
summary, this signal and system adaptive process requires two models 
and two adaptive loops. The two models are of different nature. The 
first is an ideal version of the process, and the second is of the 
same order as the process. The second model is optimized for the 
input by the first adaptive loop operating upon the difference 
between the output of the two models. The second adaptive loop 
causes system parameters to follow the optimized second model par- 
ameters, and at the same time corrects for disturbances. 

Model Output Modulation 

Dorf and Byers have shown the feasibility of the model reference 
self adaptive system using output modulation techniques. The analysis 
of the techniques and a discussion of experimental results are found 
in Ref. 5. They are briefly summarized here. Fig. 3 presents the 
block diagram of the system. 

Fig- 3 

MODEL MODULATION ADAPTIVE LOOP 




= 10 - 



In this adaptive loop, the output error is: 



E 



(t) 



(t) (t) 



1 1 + A sin tu 




2.1 



For this example error squared has been used as the criteriono Thus, 
after being shaped and demodulated a phase-sensitive doC, signal is 
obtained. When this signal is subjected to integration and filtering, 
the open loop output of the adaptive circuit becomes: 



2 

A C 



m 



(C 



m 






2.2 



The adaptive circuit loop was closed by this signal being inserted 
into the system as a driving function, to correct a system parameter. 
The result of this multiloop circuit is a com.plex, non“linear system 
which has adaptive characteristics to maintain an approximately in^* 
variant system transfer function. In functional notation the system 
output might be described as a function of the following variables: 

•^(t) “ ^ ‘“i’ hl> '"f’ 

K^d * adaptive loop gain 

= frequency of input signal 

R « magnitude of input signal 

= filter network time constants 
T * system transfer function, uncorrected 

M ~ model transfer function 



The adaptive loop filter was designed to filter out the normal system 
operating signals. The speed of adaptation depends partially on the 
filter transfer function. The gain and lag of the filter influence 
the magnitude of the adaptive response and the input signal contami- 
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nation.. One important feature of this method is that an absolutely 
constant frequency oscillator is not required. The operation of the 



adaptive loop is not dependent upon frequency stability of the modu^^ 
lating frequency. Furthermore, as there are no significant time 
constants between modulation and demodulation, the frequency of 
modulation is limited only in bounds at the lower frequency by the 
necessity of being an order above the input signal frequency. Dorf 
and Byers further investigated a two parameter correction situation, 
driving one adaptive loop with output and the second adaptive loop 
with the derivative of output. In their example the system open 
loop transfer function was: 



The model open loop transfer function was similar. The first adaptive 



loop damping (2 £ 0)^ + A 2 Q authors concluded that a two 

parameter self adaptive system of this kind was feasible, A major 
restriction was found to exist if the two parameters being adapted 
were not independent, and the adaptation would be successful only if 
the range of required adaptation was restricted. 

Application of Model Modulation 

To verify the results of model modulation, and to obtain exper- 
ience with the characteristics of the method, the method was simulated 
on a Donner 3100 analog computer. Initial investigation concerned 





2.3 



loop was driven by and to correct (K^ + A K) , the open loop gain. 

The second adaptive loop was driven by C and C to control the inner 
^ ^ m s 
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the use of a second order system with one variable parameter being 
adapted to follow a second order model. The model and the system 
were ordered for an open loop transfer function: 



Gcs'i = 






loo 



S S 1 2. S ( S I 

The parameter drift then is represented as 



2.4 






“Ur* 



w 



s(s + 2^-03^ + ) 

where n was a step parameter drift. The input was also of various 
wave forms of frequency less than the bandwidth. The adaptive system 
open loop is represented as: 



2.5 



G [o -- 



4 7i - -m^) 

where m^ is the correction derived from the output driven adaptive 
loop. This system might be considered to simulate a second order 
position servo as displayed in Fig. 4. 

Fig. 4 

SECOND ORDER SERVOMECHANISM 



2.6 




Constant Armature Current D. C. Motor Drive 
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Thus the variable drift simulated a system with varying or unknown 
triction(f), a system with constant gain but varying velocity error 
coefficient o The correction would be inserted by varying tachometer 
gain feedback magnitude. 

For operational simplicity the error criterion was altered from 
the integral of error squared, ISE, to the integral of absolute 
error or lAE. Additionally, the modulation and demodulation was 
accomplished through use of electronic multiplierso Incidental noise 
problems with these multipliers was such that the pure integration 
followed by a filter was removed from the adaptive loop, and instead 
a double filter network substituted. The alteration is shown symbol- 
ically in Fig. 5. 



Fig. 5 

MODIFIED MODEL MODULATION LOOP 




Satisfactory results were achieved with this circuit with the filter 
break points a * .1, and b = ol. The full computer diagram is shown 
in Fig. 6. The value of the gain of the adaptive loop for the circuit 
of Fig. 6 was 3- Fig. 7 presents the output data for this same circuit, 
illustrating adaptation for various parameter drift conditions 
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Additional investigation was made of the control of a third order 
system with one variable parameter by a second order models The third 
order system investigated open loop transfer function was 



G(s) 



K 100 

, ,, 

s (s+K) (s + lO + n'-m) 



2.7 



The model utilized was that of the previous example. The analog 

computer setup was that of the transfer function format similar to 

K 

Fig. 6, with an additional function p installed. The adaptive 

s + K 

P 

loop was unchanged. This system then simulated a third order control 
system as illustrated in Figo 8. 
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Fig. 8 

THIRD ORDER SYSTEM 




For this situation the variable parameter simulated would be p, 
with the correction being inserted through a small servo varying m^. 

The adaptation was satisfactory for = 40 , when p = 10. Adapt- 
ability became difficult for < 40 for large parameter shift. Fig. 9 
presents the output data for this circuit for a representative = 40 , 
under various parameter drift conditions. 

As a final illustration of the model modulation method, cancel- 
lation techniques were used for system control. The process is best 
illustrated in Fig. 10. 
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Fig. 10 

CONTROL SYSTEM CASCADE COMPENSATION ADAPTATION 




This technique is most useful when the operating system pole is not 
amenable to correction because of the physical characteristics of the 
system. In the operation illustrated by Fig. 10 zero drift caused 
compensator cancellation. When drift was introduced the correction 
forced the zero of the compensator to cancel the system pole, and the 
desired system operation was restored. The analog computer operating 
schematic is presented in Fig. 11. For operational simplicity the 
drift was inserted as a step change in the system pole. The computer 
output data is presented in Fig. 12, where the successful adaptation 
of the system is well illustrated. 
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Fig. 11 

ANALOG COMPUTER SIMULATION OF COMPENSATION ADAPTATION 




The results of these simulations of the adaptive system were 
considered to satisfactorily demonstrate the feasibility of the 
method and to serve for participant training. The analytical proof 
of the correctness of the alterations to the initial formulation, 
lAE for ISE as criterion, and filter instead of pure integration, 
have been neglected. The demonstrated accuracy of the adaptation 
was felt to be sufficient proof of the method. 
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chapter III 



IDENTIFICATION THROUGH MODEL MODULA'IION 

The methods of successful adaptive control through model roodu-- 
lation were directly applicable to the identification of operating 
systems. The requirements initially set were that the identification 
should proceed utilizing only the operating system input and output 
without disturbing the system operation. To utilize the model modu-^ 
lation method, the model would be placed in parallel operation with 
the system. The adaptive loop output would be applied to the model 
parameters to cause the model to adapt to the system. Thus, knowledge 
of the unadapted model governing equations combined with the magnitude 
of the adaptive loop output would serve to define the system. 

This method is shown in block diagram form in Fig. 13 

Fig. 13 

MODEL MODULATION WITH MODEL ADAPTATION 




This method of system identification does not relieve the engineer 
of the necessity of system analysis- The form of the transfer function, 
the order of the operating equation, and the general magnitude of 
the poles and zeros of the system must be known in order to synthesize 
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the model Ideally, one or two critical parameters whose magnitude 
may not be susceptible to analysis may be estimated, then measured 
or identified by the model modulation. These limitations should be 
expected from any adaptive system. The knowledge required to build 
an adaptive loop which would assume no knowledge of the system, and 
then synthesize perfect identification, with sufficient economy to 
be practical, is beyond the state of the art 

The standard second order servo or feedback system equation is 
formulated in transfer function terms as- 



Cal _ 

■ s’ t t 



For this system the standard open loop transfer function is; 



3. 1 



3 . 

The response of the second order system is completely described by 

the two variables, zeta and ou . A higher order system, or any 

n 

system with closed loop transfer function zeros, naturally, has 
additional significant variables. If, however, it is considered that 
there are two variables of primary interest, the higher order open loop 
transfer function may be written as: 

3 , 

^ s(s+2i-w;) # ("sip ) 



where cu and zeta are these variables, 
n 
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In the development of model modulation system ident i f i ca tion method 
two major assumptions were made. The tlrst of these Is. 

1 

That in the operating system to be identified there will be but 
two unknown parameters. Thus the system must be second order, 
a completely dominant second order approximation o£ a higher 
order system, or a higher order system which has the ether 
parameters fully identified If the latter is the case, the 
format of the preceding equation will be usted to df^^:cribe the 
system. 

This first assumption limited the parameter adaptation problem to the 
two variables zeta and 



Examination of any physical system reveals that zeta and may 

or may not be functionally independent The simple second order system- 
detailed in Fig. 5 serves as an adequate example. The closed loop 
transfer function is: 



In this system variation of the friction(f) will yield an independent 



n 




3.4 



zeta and (ju „ Variation of all other factors will result in zeta and 



n 



(D being different functions of the same variable, and thus not in- 
n 



dependent. Thus zeta and o) may be independent or both different 



functions of the same variable. 
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The adaptive method that has been discussed to this point has 

dealt with but one variable. Since the interdependence of zeta and o) 

n 

is at the most one of their being two different functions of one factor, 
it is reasonable to expect that two adaptive loops will be required 

The second major assumption required depends heuri s tica lly on this 

discussion of the functional relationship of zeta and ud , This 

n 

assumption is; 

IIo 

That model may be adapted to the system through two separate 
adaptive loops. The first of the loops will operate upon the 
system and model output to develop the first corrective factor. 

The second adaptive loop will operate upon the derivative with 
respect to time of the system and model output to develop the 
second corrective factor. 

The remainder of this report is devoted basically to showing that these 
two assumptions are justified 



Proceeding on these assumptions the investigation utilized models 
of the following general form: 



Gcs.^ 






3.5 



A schematic of the operation of these loops as tentatively applied to 
system identification is illustrated in Fig. 14. 
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Fig. 14 

SCHEMATIC OF TWO PARAMETER CORRECTION 




With this model the characteristic equation of the operating 
system could then be easily detailed. For the model with the form 
of equation 3.5 and unity feedback, the closed loop transfer function 
becomes : 




This function may then be used to accurately detail the system differ*^ 
ential equations. 



Only limited single parameter, single adaptive loop correction 
investigation was made. Primary investigation concerned the two par 
ameters of the previous equation as in Fig. 14. 

Figc 14 also illustrates the visualized mechanization of the model 
modulation identification system. Potentially the model and the adaptive 



- 32 ^ 



loops would be constructed through use of standard operational amplifier 
networks. Thus any moderate sized analog computer operating in con- 
junction with the scaled input and output of a physical system would 
serve as an identifier. 
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Chapter IV 



ANALOG COMPUTER INVESTIGATION 



Model Modulation Identification o l an Equation to Order 

The primary investigation was centered around use of the analog 
computero The investigation concerned use of two undetermined par- 
ameters in the system. The simulated control system transfer function 
was of the form shown In formula 3.3 The model transfer function as 
of the form shown in formula 34. In both cases the actual computer 
arrangement was that of the transfer function format. The system 
values of and zeta were set up through switching circuits so that 
three different values of each could be used directly The adaptive 
loops used on the computer were those shown in Fig 6 In general, 
both filters in the adaptive loops were set for the break at 
omega = .1 radian. Fig. 15 presents the full analog computer diagram 
for a simulated third order system with a model tracking two unknown 
parameterSc This diagram conforms to the format of formula 33 and 35. 
The initial investigation with this computer arrangement concerned a 
control system with simulated open loop transfer function of: 



G(s) - 



10 (100) K 
s (s -f p) (s -f 100) 



5 ^ K ^ 15 
5 ^ p ^ 15 



4. 1 



and a model transfer function: 



G(s) 



1000 (10 m ) 

' 

s (s + 100) (s 10 +m») 



4.2 



34 
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Analog' C!5aiputer DiPfrao) for Two Parameter 
Mod el Ada p 1 8 1 1 on 



The input was primarily deterministic. It consisted of a square 



wave with a frequency which was less than the bandwidth. Tests were 

also carried out however, using a random input. This random input 

was the output of a Gaussian **white" noise generator with an 18 db 

cut off at approximately the system bandwidth. At the completion of 

each run, the recorded m^ and m^ were used in conjunction with the 

known model characteristics to derive the closed loop characteristic 

equation of the simulated system. This allowed comparison of actual 

system zeta and U) with those derived from that identified as the 
n 

system. A non-dimensional corrective parameter was examined and 
used as a partial criterion of identification. This corrective 
parameter a^ is best understood through examination of Fig. 16. 



Fig. 16 

CORRECTIVE PARAMETER 
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In this figure, is the dominant root location of the simulated 

control system, (2) is the root location of the uncorrected model, 
and is the position of the corrected model root. The arithmetic 

mean ¥ and the standard deviation, or of were computed for each 
series of tests. 
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The data for the first series of tests is shown in Table 1. In 



this table the variance of the tracked parameters is shown. The 

roots, zeta and U) have been shown for only the two extreme cases 
n 

in each test, that of the maximum and minimum zeta. The correction 
parameter a^ and the sigma have been calculated from the assembled data 
of the eight situations for each test where the parameters do not 
match. Fig. 17 shows a sample of the analog computer output. This 
data is that of the first and last cases tabulated in Table 1 
Figure 17D will be used to illustrate the method of computing the 
characteristic equation roots of the operating control system. For 
this example, the model closed loop transfer function was 

^ 10.7 (s+14) (10*m^) 

R " s (s+15) (s+104mp + 10.7 (s+14) (10-m^) ^ 

From figure 17D the values of the correction factors may be seen to 

be m = 1.5 and m* = 1.080. The characteristic equation is then 
c c 

+ 26.08 + 289-1 s+1724. 4 4 

The roots of this equation were obtained through digital computer 
service. These roots are: 

-5.894 + j 9.268 

-14.29 

Fig. 18 shows a plot of the characteristic equation dominant root 
positions for the actual system and the identified system detailed 
in Fig. 17. 
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Equation Order Approximated 



During the process of analyzing an operating system, it is 
inevitable that approximations would be made. These approximations 
will generally mask or ignore the very small time constants, which 
in the s plane are the poles of large magnitude. In trying to use 
the model modulation method to identify an operating system it is 
then inevitable that the model and the system would be of different 
order. For a step input, the effect of these poles would be greatest 
during the initial portion of the transient response, and the effect 
would become inconsequential in the later portions of the response. 

It was hypothesized that if the uncertain parameters of the operating 
system concerned dominant roots, satisfactory identification would 
take place. The restriction to this hypothesis was two fold. It 
should be expected that if the input was of a random noise nature 
with no stead state, no identification would be possible unless, as 
in the first test of Table I, the equation order and the additional 
parameters are known. Additionally, with deterministic inputs, 
identification could take place only when these unknown poles were 
of sufficient magnitude to be unimportant with respect to the dominant 
response . 

To examine this restriction the tests detailed in Table II were 
hold. Since the systems to be tested were set up in transfer function 
format on the analog computer, each pole was a new block set into the 
basic arrangement of Fig. 14. The information in Table II is of the 
same nature as in Table I, but as all nine system configurations differed 
from those of the model, all nine situations were used to calculate 
the a and o. 
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LOCATION OP DOMINANT ROOTS K)R THE 
ACTUAL AND IDajTIFIED SYSTEMS 
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The results of these tests were as expected The random input 
revealed, as expected, lack of identification. Under deterministic 
inputs, when the postulated restrictions were held, satisfactory 
identification took place. One interesting point revealed in this 
investigation was that an adequate test for satisfactory identification 
would be made if the frequency of the input could be varied. Naturally 
this violates the basic restriction that led to the development of 
this method, i.e., no violation of the natural operation of the system, 
but is a point of interest. For this process the frequency was varied 
between a lower frequency with a period equal to the second order 
time constant and an upper frequency equal to the bandwidth. When 
the adaptive loop outputs remained relatively constant through this 
frequency range it was found that these outputs served to adequately 
identify the dominant roots when the characteristic roots were founds 
If, however, the identification varied through this range, the effect 
of the ignored poles precluded satisfactory identification. 

Figure 19 presents illustrative data from the first tablulated 
case of Table II for this series of tests. Fig. 20 shows the dominant 
roots of the system for each of the tests illustrated in Figc 19, 
with the dominant roots as identified from the data. 

Examination of the data reveals that for both of these types of 
equation situations good dominant root identification has been achieved 
It was found that these dominant roots could be identified to within 
an average error of 2 . 27 o of the undamped natural frequency and to 
within an average error of 4.76% of the complex phase angle. When 
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this method is used to identify the characteristics of an actual linear 
control system, it is highly probable that greater accuracy may be 
attained since the model parameters may be varied to obtain error 
bracketing. This process is illustrated in Chapter VI 
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Chapter V 



DIGITAL COMPUTER SIMULATION 

Direct utilization of a digital computer for identification of 
an operating control system was not considered 1 t is possible, 
however, that operating systems might be monitored or controlled by 
digital techniques. In this case the output would be presented in 
the form of discrete numerical data. Ideally this data could be used 
as input to a digita 1- to-analog converter, and the analog techniques 
previously described could be used for operating system identification. 

In the absence of this equipment, digital simulation of the control 
system, model, and adaptive loops was investigated to illustrate the 
techniques and the potential applicability of model modulation techniques 
through numerical techniques. This investigation was conducted through 
the use of the CDC 1604 Digital Computer. The program code was written 
in Fortran. A standard Runge-Kutta numerical integration process was 
used to solve the differential equation which governed the operation 
of the full model, system and adaptive loop. The differential equations 
were formulated in a transfer function format. 

Three basic configurations of model adaptation were tested through 
this digital computer simulation. For these configurations the model 
was used to identify the operating control system characteristic 
equation roots under the following conditions 

1. Second order control system with undetermined pole. 

2. Second order model, third order control system with uncertain 
gain and uncertain dominant pole. 

3. Third order system and model with two uncertain poles. 
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The actual investigation was limited to achieving a successful solution 
for one parameter set for each of these configurations. The excessive 
time used by integration techniques on the digital computer was of 
course a limiting factor. The first of these configurations will be 
by-passed in this discussion since it can be considered to be a sub- 
system of the second configuration. 



The model for the second configuration had the following open 
loop transfer function: 



10 (10 - m ) 
c 

" s (s + 10 + m.) 

c 

and a closed loop transfer function of 

^ 10 (10 - m ) 

^ V „ ^ 

R '®''m ~ s (s + 10 + m.) + 10 (10 - m ) 

c c 



5.1 



5.2 



The control system open loop transfer functions were 



G (s'^ 



IS.OOO 

S (“sf loo)(s+8) 



5. 3 



and 



(s') ^ 5 4 

S ( s^- ioo)(z4~ 

Figure 21 presents the transfer function format which was used to 
formulate the differential equation for the FORTRAN program. This 
format allowed formulation of the operation by a set of quasi-state 
variables of the form: 
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The solution of the equations was programmed to be presented in both 



numerical and graphical format. Fig. 22 presents a sample of the 
input program Fig. 23 presents a sample of the numerical data output, 
and Fig. 24 presents the graphical data output for the control system 
open loop transfer function given by equation 5 3. 

The graph of Fig. 24 presents the corrective parameters as 
"corrected pole" = (10 + m.) and "corrected gain" = (10 - m^) 

Table III shows the identification of the characteristic equation 
roots for these two situations. 



Table IIIA 



DIGITAL COMPUTER CHARACTERISTIC EQUATION ROOTS 



Control System 



Model as Corrected 



Equation 5. 3 



3.36 + j 10.35 
101.27 



3.07 + j 9-97 
101.15 



Equation 5.4 



5.55 + j 6.96 
100.9 



4.88 + j 7.77 
100.93 



- 62 - 



I 



Fljurp 2.2. a. 



Dii'ltRl lonputfir Program 



C 

c 

c 

c 



PROGRAM ADAPT 12 

N=NO.OF EUNS., NP=NO. OF INCREMENTS BETWEEN PRINT OUTS, TO=INITIAL 
TIME, TF=FINAL TIME, CT=TIME INCREMENT, TN=TIME FIRST PRINTED DATA, 
X(J)*INITIAL VALUES OF THE VARIABLES X(J) 

NG=NO. OF INCREMENTS BETWEEN GRAPH POINTS 

GDI MENS ION X ( 30) ,Xl ( 900) ,X2 (900 ) ,X3{ 900) ,X4( 900) , X5{ 900 ) , X6 { 900 ) , T I 
1ME(900), H(30) 

COMMON R, E, T2, H 
100 FORMAT (31 1C/IPF10.0) ) 

102 FORMAT { 12 / ( I2,F2C.O ) ) 

200 FORMAT (F8.2,1CF11.5) 

2010F0RMAT (99H TIME SYS.POS. MOD. PCS. POS.DIFF. SYS.VEL, 



GAIN COR. POLE COR. 



//) 



2H). / ) 

, 27H, POLE CORR. 



IMOD.VEL. VEL.DIFF. 

202 FORMAT (/) 

203 FORMAT (IHI ) 

204 FORMAT {//) 

205 FORMAT {^9H SYSTEM^ F9. 1, 1 3H/S ( S+1 00 ) ( S+ F6.2, 

2060F0RMAT (,32H GAIN CORR FILTER CORNER AT W = 0.6»2 

1 FILTER -AT W = F6.2 /) ^ 

207 FORMAT f25H MODEL = 1 CO . /S ( S+ 1 0 ) . V/ ) ’ ■ 

2080F0RMAT f57H ADAPTIVE SYSTEM GAINS, OUTPUT (S) = F6.2, 20H, FIRST DE 
IRIV. (01' = F6.2 /) 

210 FORMAT (20H COEFFICIENT NUMBER 12, 2H = F12.4, IH. ) 

READ 100, N, NP, NG , TO, TF, OT , TN, (X(J), J=1,N) 

READ 102, NF, ( K, (H(K)), I = 1,NH ) 

PRINT 204 

PRINT 208, F( 16), H( 17) 

ASYSG = H(5)* H(2)* H(3) 

PRINT 205 , ASYSG, H(4) 

, H( 1 1 ) , H( 15 ) 



I = 1,NH 



1 



PRINT 206 
PRINT 207 
DO 6 

PRINT 202 
PRINT 210, 
PRINT 203 
PRINT 201 
T = TO 

T2 = 0.0 
NUMPTS = 0 
MPTS= 0 
DPOS = Xtl ) 
OVEL = X(2) 



I , H ( I ) 



X(U) 

X(5) 

AGACOR = H(7) - X(7)*H( 16) 

APCLCOR = H(6)*( 1. + X(9)*H(17)) 

IF (T-TN) 36, 2, 2 
IF(MPTS) 3C,35,30 
IF ( XMOOF (MPTS, 50*NP ) ) 31,33,31 
IF ( XMODF (MPTS, 10»NP) ) 32,34,32 
IF( XMODF( MPTS, NP ) ) 36,35,36 

33 PRINT 203 
PRINT 201 

34 PRINT 202 

35 PRINT 200, T , X ( 1 ) , X ( 4 ) , DPOS , X ( 2 ) , X ( 5 ) , DVEL, AGACOR , APOLCOR 

36 MPTS=MPTS+1 

37 IF ( XMODF (MPTS, NG ) ) 38, 39, 38 
39 NUMPTS = NUMPTS+1 

" 40 



2 

30 

31 

32 



IF (900 - NUMPTS ) 38, 38, 
40 XI ( NUMPTS) = X ( 1 ) 



X2(NUMPTS) 

X3( NUMPTS ) 
X4(NUMPTS ) 

X5( NUMPTS ) 

X6( NUMPTS) 
TIME (NUMPTS) 



X(4 ) 
AGACOR 
APOLCOR 
DPOS 
DVEL 
= T 



38 IF(TF-T-DT) 10,20,20 



10 PRINT 204 



CALL 


GRAPH 


( 


NUMPTS 


f 


T IME 


9 


X3 


9 


8 


) 


CALL 


GRAPH 


( 


NUMPTS 


f 


T IME 


9 


X4 


9 


8 


) 


CALL 


GRAPH 


( 


NUMPTS 


9 


T IME 


9 


X5 


9 


8 


) 


CALL 


graph 


( 


NUMPTS 


9 


T IME 


9 


X6 


9 


8 


) 


CALL 


GRAPH 


( 


NUMPTS 


9 


T IME 


9 


XI 


9 


8 


) 


CALL 


GRAPH 


( 


NUMPTS 


9 


T IME 


9 


X2 


9 


8 


) 



STOP 

20 CALL RKUTTA(N, T, X, DT ) 
T2= T2 + DT 
T = T + DT 
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Fisru rfi 2.2. h 



Digital 'Computer Program 
( con ' t ) 



C 

C 



C(»4) , H(30) 



tNU 

SURRCUTINE RKUTTA (N,T,X,0T) 

DIMENSION X(30) ,AK(U,30) , XDCT(30), XC(30)t 
COMMON R, £, T2, H 
C ( 1) =0.0 
C(2) =0.5 
C( 3)=0.5 
C( 4 ) = 1 .0 
DO 4 1=1,4 

TC = T + d ( I ) « OT 
DO 2 J=T,N 

XC(J)= X( J) +cm* AK(I-1,J ) 

CALL DERIVC TC,XC,XD0T ) 

CO 4 Jf 1 ,N 

AK( I ,J)= DT« XCOT(J) 

DO 3 J=1,N 

X( J)=X(J)+( AK(1 ,J)+2.»AK(2,J)+2.* AK ( 3 , J ) +AK ( 4 , J ) ) / 6. 

RETURN 

END 

SUBROUTINE CERIVIT, X, XDCT) 

DIMENSION XDOT(30), X(30), H(30) 

COMMON R, E, T2, H 
IF( 1 .0- T2)3,3, 1 
R = H( 1 ) 

GO TO 2 
R = -H(1 ) 

IF(2.0 -T2) 4,4,2 
T2= 0.0 
GO TO 1 

E = SINF(314.*T) 

H(4) = VARIABLE PCLE 
H(5) = VARIABLE GAIN 



XDOT( 1 ) 

XDCT (2 ) 

XD0T(3) 

XD0T(4 ) 

XD0T( 5 ) 

X00T(6) 

XD0T(7) 

XD0T(8) 

XD0T(9) 

RETURN 

END 

END 



H(5)«X(2) 

X (3 )«H( 3) - X (2) «H(4 ) 
H(2)*(R - XM )- X(3) ) 

- X(7)*H( 16)) 



n \ ^ 9 ^ \ r\ — y % f — \ ^ f f 

X(5)«(H(7) - X(7)*H( 16)) 

= H(6)«(R - X(5) - X(5)*X(9)*H( 17) - X(4)) 

= H (8 )*E*A8SF(X ( 1 )- X(4) - E*X(4)) - X(6)*H(9) 

= H(10)»X(6) - X(7)«H( 1 1) 

= H ( 12 )*E*ABSF( X( 2) - X(5) - X(5)*E) - X(8)»H(13) 
= H(14)*X(8) - H( 15)*X(9) 
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Figure 2-3 a 

Dlffltal 'l'):nputf?r Nuuwrlcal Output 



ADAPTIVE SYSTEM GAINS, OUTPUT (S» = 10. 00, FIRST OERIV. (0) = 

SYSTEM= 12000. 0/S(S+lC0) ( S+ 8.00). 

GAIN CORR FILTER CORNER AT W = .10, POLE CORR. FILTER AT W 



MODEL = 100. 


/S I S+ 10 ) . 
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NUMBER 
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i 



The third configuration examined was that for the following model 



open loop transfer function: 



G(s) 



10,000 

s (s + 100 +m ) (s + 10 +m.) 
\ c' ' c 



5 6 



with control system open loop transfer functions of; 



G(s) 



10,000 

s (s + 150) (s + 15) 



5.7 



and 



G(s) 



10,000 

s (s + 75) (s + 7.5) 



5.8 



The equations were formulated in the manner of Fig. 20, except 
that the adaptive loops each utilized one filter and one integrator 
instead of the two filters previously detailed. Satisfactory results 
were obtained. Fig. 25 presents sample results from the investigation 
of equation 5.7, and Table IV presents the root identification 
achieved. 



Table IIIB 

DIGITAL COMPUTATION CHARACTERISTIC EQUATION ROOTS 
TWO UNCERTAIN POLES 



Equation 5. 7 



Control System 

- 7.25 + j 3. 72 

- 150.49 



Model as Corrected 

-7.87 + j 2.1 
-150.49 



Equation 5.8 



- 2.81 + j 11.05 

- 76.87 



-3.45 + j 11.05 
- 74.56 
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Digital Goaiputer Graphical Output for Integration in Ada.ptive Loops 



Chapter VI 



ACTUAL SERVOMECHANISM IDENTIFICATION 



Intro du e t ion 

The previous discussion and the data in Table I and Table IT 
indicate that the identified location of the roots of an operating 
system is the function of many variables. The type of model and 
the proximity of the model equation roots to the roots of the equation 
of the operating system are obviously factors of importance in this 
identification. With this point in mind, it is suggested that the 
following steps be followed in the process of identifying an actual 
system. 

1. Analyze system operating factors to enable detemina tion of 
open and closed loop transfer functions. 

2. Determine which values in the open loop transfer function 
are questionable. 

3. Utilize operational amplifiers or an analog computer to 
prepare the dynamic model and the adaptive loops. 

4 Scale the model input and outputs to the estimated 
corresponding values of the system. 

5. Set indicated adaptive loop polarities to adapt the model 
to the undetermined or questionable values in the system. 

6. Adjust gain of adaptive loop to obtain satisfactory operation 

7. Calculate the transfer function of loop corrected model, 
scale a new model to these functions and repeat operation 
to generate better conformance or to confirm values pre- 
viously attained 
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A position servomechanism was used to illustrate the applicability 
of the model modulation techniques to an actual operating system. 

The servo used was one developed for use in instructional laboratories 
by Professor J. R. Ward of the U. S, Naval Postgraduate School, and 
is partially described in Ref. 6. The full description and details 
of the operation of this system are on file at this school. 

Figo 26 presents a block diagram of the system. This repre» 
sentation shows that the servomechanism may be considered as a gear- 
head system. Thus the use of this identification method conforms to 
the stated requirement of parallel model operation. 

Fig. 26 

POSITION SERVOMECHANISM 




Analysis of this mechanism resulted in the transfer function diagram 
of Fig 27. 
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Fig. 27 

POSITION SERVOMECHANISM BLOCK DIAGRAM 




The approximate parameter values in the system were determined. The 
magnitude of inductance (L) was such that the approximation L 0 
was made. The open loop transfer function was then; 






G(si =■ 



RJ" 



•m 






S 4- 



RJ / 



)< 

sCsi- p) 



6 1 



The value of friction(f) was difficult to evaluate because of the 
expected non-linearity of friction in a d.Co motor. A nominal value 
was chosen to represent the viscous friction, while coulomb friction 
effects were ignored. It was further assumed that the open loop gain 
was questionable. As a first estimate it was assumed that the open 
loop transfer function was: 



r - k ^ 400 

' s (s4 fT sCs+8) 

Fig. 28 presents the control sys tem-model -adaptive loop block 
diagram. 
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Fig. 28 

SERVOMECHANISM IDENTIFICATION BLOCK DIAGRAM 




The actual model and adaptive loop analog computer schematic was that 
detailed in Fig. 15, while the scaling factors were achieved through 
cascaded operational amplifiers to assure correct polarity. The 
servomechanism and the model were driven by a square wave of one cycle 
per second. Table IV presents the results of four tests made first 
to bracket, then to closely approximate the open loop transfer 
function. Fig. 29 presents the results of the final run listed in 
Table IV. One the basis of these results the transfer function of 
this servomechanism was estimated to be: 



+ lO.IS's + G2i? 

These results compare favorably with information derived from other 
sources. Calculated results derived for careful component analysis 



6 . 3 
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Table IV 



POSITION SERVOMECHANISM CHARACTERISTIC EQUATION IDENTIFICATION 
OPEN LOOP TRANSFER FUNCTION SYSTEM DYNAMIC EQUATION 



Model 


Corrected Model 


Identified Roots 


10 (40-m ) 


574 


-4.83 +j 23.45 


c 

s (s+8+m^) 


s (s+9.66) 


C = .201 OD =24 
n 


10 (70-m ) 
c 

s (s+10“hn*) 
c 


652 

s (s+10.15) 


-5.075 +j 25 

C = .204 UJ =2 5,55 
n 


10 (55-m ) 
c 

s (s+lO+m^) 


622 

s (s+10.24) 


-5.12 +j 24.2 

C = .206 UJ = 24.95 
n 


10 (60-m^) 
s (s+10"hn^) 


625 

s (s+10.15) 


-5.075 +j 24.4 

C = .208 UJ =25 
n 
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indicated a closed loop transfer function of 



fs) 



S9o 

s* + (^ n.s t i.s)^ 4 59o 



6.4 



while various open loop frequency response tests indicated the closed 
loop response to be: 



. 3-6^ ±70 

4 (lot i)% 4 s(>s-± no 

In view of the non-linearities involved in a d.c. servo motor it was 
felt that these model identified roots represented a successful 
illustration of the validity of the method. 



Discussion 

These servomechanism tests confirmed that the model modulation 
method could be successfully applied to obtain the governing operating 
equations of a servomechanism. Beyond this the tests emphasized the 
necessity of continuous excitation of the control system for successful 
model adaptation. For a linear control system the model modulation 
adaptation will operate only while the servo is undergoing excitated 
or active response to an input signal. It would be expected that the 
correction or adaptive signal would be generated in response to a 
step for a period equal to the settling time of the system. Beyond 
the settling time, with no excitation response operation to generate 
error, the time constants of the adaptive filter network will allow 
dissipation of the adaptive signal magnitude. The analog computer 
simulations detailed in this chapter illustrate this operation of the 
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adaptive loop For the analog tests the general input frequency of 

the square wave was 0.5 cps. The average settling time for the 

second order model was T — 4/6u3 1.0 sec, and the bandwidth was 

n 

f 1.9 cps When tests were run with systems with small zeta , lower 
input frequencies would achieve successful adaptation. In all cases, 
successful operation was attained with frequencies which were as high 
as the bandwidth frequency. 

In the case of the operating servo, the non-^linearities of the 
friction caused a faster damping than the linear interpolation would 
indicate. For this system, the adaptation proceeded successfully 
when the driving square wave frequency allowed limited oscillation 
per step. This in effect was a partial masking of the non-linearity 
through testing techniques. 

These situations emphasize the limitation of control system 
excitation in the use of the model modulation method It is obvious 
that this limitation exists for both system adaptation and for system 
identification through model adaptation A potential solution of 
this problem might be found through investigation into utilization 
of signal identification circuits. This signal identification might 
be applied to the output of the control system to control the operation 
of the adaptive circuit. This identifier could allow the adaptive 
circuit to operate while the control system would be underdoing 
oscillatory or excited response. When this control system response 
became quiescent, the identifier could place the adaptive loop into 
a hold state with a constant level of output maintained to maintain 
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the adaptation until such time as the control system would be again 
activated. Additionally, it is expected that the adaptive loop 
elements would be synthesized through the use of operational ampli** 
fiers. Since these are magnitude sensitive, another role of the signal 
identifier would be to hold adaptive loop operation when signal mag- 
nitude became excessive. 

The techniques of sampled-data control systems offer attractive 
possibilities for circuits to fulfill these requirements. Hybrid 
digital-analog circuits might also be investigated. In any case, 
these techniques warrant additional investigation. 
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Chapter VII 



CONCLUSIONS AND RECOMMENDATIONS 



Cone lusions 

It is concluded from the results of this investigation that the 
identification of the characteristic equation roots of a control 
system by model adaptation through model modulation is a valid 
technique. The technique is subject to operating restrictions as 
follows : 

1. The control system must be linear or be subjected to operating 
non-linearities of minor magnitude. 

2. The open loop transfer function must have but one or two 
parameters whose magnitude is uncertain. 

3. The additional parameters of the open loop transfer function 
must be known or be of sufficient magnitude to be incon- 
sequential. 

4. The control system must be undergoing active operation in 
response to an input signal. 

It is further concluded that: 

1. The roots of the characteristic equation of a control system 
may be identified with a high degree of accuracy. During 
these tests, these dominant roots were identified to within 
an average error of 2.27o of the undamped natural frequency 
and to within an average error of 4.767o of the complex phase 
angle . 
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2. When the order of the characteristic equation is unknown, 
second order dominant mode approximations will yield 
consistent results which are valid in the approximation 

Recommendations 

The results of this investigation indicate that there are 
certain areas of interest which warrant further investigation. The 
first of these is investigation into the stability of the dynamic 
adaptive loop. The stability limits are subject to various major 
bounds in loop gain, input signal frequency, and input signal magnitude. 
Theoretical analysis of these bounds will be extremely difficult, 
thus analog computer investigation might offer a valid method of 
investigation. 

The second of these would be the investigation of methods to 
remove the necessity of a continual excited response from the control 
system. Certain suggestions have been made in the text which might 
be applied to achieve this end. 

The third suggested area of investigation is in the identifi* 
cation of control system operation subject to known non-linearities 
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